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Abstract 

A "Kallen-Lehman" approach to Ising model, inspired by quantum 
field theory a la Regge, is proposed. The analogy with the Kallen-Lehman 
representation leads to a formula for the free-energy of the 3D model with 
few free parameters which could be matched with the numerical data. 
The possible application of this scheme to the spin glass case is shortly 
discussed. 
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Introduction 

The three dimensional Ising model (henceforth 3DI) is one of the main open 
problems in field theory and statistical mechanics. A large number of interest- 
ing statistical systems near the transition point are described by 3DI and the 
theoretical methods suitable to deal with such a problem manifest deep connec- 
tions in various areas of physics ranging from quantum information theory to 
string theory (to provide with even a partial list of references is really a com- 
pletely hopeless task, a detailed review- which, by the way, does not cover all the 
approaches and results on the subject-is [IT]). Besides its intrinsic interest is 
statistical physics, since the formulation of the Svetitsky-Yaffe conjecture [T3] , 
it has been widely recognized its role in describing the deconfinement transition 
in QCD. For this reason, the 3D Ising model is worth to be further investigated. 

Here, a " phenomenological'f]] approach is proposed which allows to formally 
connect the 3DI model to the Ising model in one dimension less. The method 
is inspired by quantum field theory a la Regge and by the Kallen-Lehman rep- 
resentation. Recently (see [2]) this approach has been able to shed some light 
on the physical origin of non extensivity. 

The paper is organized as follows: in the second section the "Kallen-Lehman 
interpretation" of the 2D Ising model will be discussed. In the third section the 
3D model will be described. In the fourth section the high temperature behavior 



The meaning of the term " phenomenological" in this context will be more clear later on. 
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will be analyzed and it will be shortly pointed out how to apply the formalism 
to the spin glass case. Eventually, some conclusion will be drawn. 



1 Ising Kinetic term(s) 



In this section the main idea of the method are introduced. 

The partition of functions of Ising models in various dimensions all belong 
to the following class of partition functions 



Z(Jab) 



E 



exp 



_ ab 



±1 Va, a,b=l,..,N 



(1) 
(2) 



where the sum is over all the spin configurations, N is the total number of spins 
in the system and the "kinetic matrix" J a b identifies the models. A useful way 
to write the Ising kinetic matrix is as follows. Let us begin with the 3DI model 
on a finite cube of n x n x n size (where n> 1): one can enumerate the spins 
starting from the leftmost spin on the lowermost face of the cube and continuing 
from left to right in any row, from the lowest row to the highest row in any face, 
and from the lowest face to the highest face in the cube. In such a way, the 
nearest neighbors of the spin number k are the spins with numbers k ± 1, k ± n 
and k ± n 2 . So that the kinetic matrix of the 3DI model can be written as 
follows: 



J, 



[f = (3 [S(\a - b\ - 1) + S(\a - b\ - n) + S(\a - b\ - n 2 )] , 
(3 = e/K B T, e>0 



(3) 
(4) 



Kb being the Boltzmann constant and e > in order to have a ferromagnetic 
interaction. From this construction it is also clear that the kinetic matrices of 
the 1DI and 2DI models arc: 

7(1) 



ab 

r(2) _ 



06(\a-b\-l), 



J K ab = (3[S(\a-b\-l) + S(\a-b\-n)} 



(5) 
(6) 



At least at a formal level, one can see the kinetic matrix of a D-dimensional 
Ising model as a perturbation of the kinetic matrix of the (D-l)-dimensional 
Ising model in which a further off-diagonal term has been added. This point of 
view is useful since already the relation between the 2DI and the 1DI models 
fits with this scheme. To see this, it is worth to write down the well known free 
energies of both models: 



F 1D {f3) = log 2 cosh (3, 
F 2D ((3) = F 1D (2(3) + ±- j dxlogl^ 



7T 

J efelogji l + \Jl- k(f3) 2 sin 2 x j 



(7) 
(8) 



HP) 



cosh 2(3 coth 2(3 



, < k{(3) < 1. 



(9) 
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The above interpretation of the 2DI model as a perturbation of the 1DI model 
is particularly clear: the free energy F^d appears as a "dressing" a la Kallen- 
Lehman (see, for example, [T5]) of the partition function Fuj. It is worth to 
recall here that the singularity in the specific heat is determined by the region 
in which k (which multiplies the sin function inside the square root) is near to 
1. 

To be more specific, in a scalar interacting quantum field theory in which 
the bare propagator is 



A^m ) = — g r> 

the dressed propagator A^ can be written as the convolution of the bare prop- 
agator with a suitable energy function p: 



A = / p^)A^(p z )dp z . 

Usually, the function p^ has S contributions (coming from single particles states) 
plus other terms (coming from the continuous spectrum) 

Pcpiv 2 ) ~ s ^ 2 ~ am ^ + Mm 2 )- 

i 

In very much the same way as it happens in field theory (in which the bare 
propagator is dressed by the interactions which lead to the exact propagator 
which can be expressed as a convolution of the bare propagator with a suitable 
energy-density function), in this case F%d appears as a convolution of F\d with 
a suitable density pgi 



F 2D (f3) = / P/3 (p)F 1D (p)dp, (10) 
Jo 

Pfl Qi) = 5(/x-2/3) + ^(/x) (11) 

7T 

= ^y#J(^|cosh/i-i [l+(l-fc(/3) 2 sin 2 0) Q ]|, (12) 



(13) 



where J{p) is the Jacobian needed to have a well-defined (5-function, the function 
k(/3) is defined in Eq. ©. Thus, the above arguments lead to think that a 
similar representation could also hold in the case of the 3DI model in terms of 
the 2DI model. 

It is worth to stress that the above dressing function ctL"^ (p) can be thought 
as an infinite sum of contributions of "moving poles": since 



S(x) ~ limP 



7 



7^0 (_ X 2 + 7 2 
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(where P {} denotes the principal valne) one can write 



7T 

<7^(/i)~ lim^P J #J( M ) 



7 2 + {cosh^-± 1 + (1 - /c(/3) 2 sin 2 </>)" j 



which looks like the sum of contributions of particles with (^-dependent mass. 
The (^-dependent "trajectories" of the (generically complex) poles are deter- 
mined implicitly by the following equations 



<] cosh // — i 



1 + (l - k{[3) 2 sii 







(where, for the 2D Ising model, a — 1/2) and one has to sum (actually, one 
has an integral over (j>) over the "moving poles". These "moving poles" are 
indeed a characteristic features of Regge theory [12] [13] . The strength of the 
Regge approach to particles physics is its generality [4] [17]: it is based upon 
very general assumption, it does not lie on perturbation theory and most of its 
results (such as the Veneziano amplitude [16 ) con be derived without a precise 
knowledge of the microscopic interaction. It allows to describe in terms of few 
parameter^] a huge number of observations. The application of Regge theory 
to statistical mechanics has been proposed for the first time in [2]. 



2 The " Kallen-Lehman" free energy 

In the previous section it has been suggested that the free energy of the 3DI 
model should be seen as a dressing of the free energy of the 2DI model. At a 
first glance, there is no reason to expect that the dressing function leading to 
F313 could bear any resemblance with the dressing function in Eq. (|11|) (which 
leads to i^u)- However, the lift from the kinetic matrix of the 1DI model to the 
kinetic matrix of the 2DI model is very similar to the lift from the kinetic matrix 
of the 2DI model to the kinetic matrix of the 3DI model. Thus, the function 
which dresses the 2DI model giving rise to the 3DI model could be quite similar 
to the function which dresses the 1DI model giving rise to the 2DI model. One 

2 Such parameters in the case of the old-fashoned S-matrix approach to the strong inter- 
actions, were deduced by experimental data [4] |17| but in principle could be computed in 
the underlying fundamental theory. The strength of the Regge approach is that it says many 
things even if one is not able to compute the above few parameters. 
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is lead to the following ansatze for the free energy of the 3DI model: 



F 



3D 



p {p)F 2D {p)dp 
A 



(14) 



7T 7T 



d<j) / dtp 



o o 













i - 


I* 







A(0) 



. 2 

sin ip 



A(0) 



(l+(l-fc(/3) 2 sin 2 0) a ) 2 , 1 < (A(0)) 2 < 4 



JAU) - 1 



(15) 



(16) 



where pj(p) is defined in Eq. (fTTj) , the function k((3) is defined in Eq. ^j. In 
a sense, from the point of view of Regge theory, the above dressing is the most 
"gentle" way to go beyond the integrability of the 2D model: one introduces 
a more complex analytic structure but in a "educated" recursive manner. The 
present field theoretical point of view only suggests an ansatze and the correct 
scheme should be to leave v and A as a free parameters to be fitted on numerical 
data. It will be shown that in order to describe the critical behavior it is enough 
the parameter a, the role of the other parameters will be also discussed. It is 
worth to note here that in this case also the factor which multiplies sin 2 <p inside 
the square root fulfils the condition (fTo| analogous to the condition © which 
gives rise to the singularity in the 2D case (see [5]). The qualitative behavior of 



F 



(i/,A,a) 



3D 



(/3) both at low and at high j3 is similar to the two dimensional case. 



2.1 The specific heat 

In the present parametrization (which is the analogous of the classic parametriza- 
tion of the two dimensional model [S]), the specific heat is proportional to the 
second derivative^ of F^^ with respect to (3: 

^3 ( D A=1 V)~C ( 2)+C (3 ) 

3 The first derivative of the free energy has no problem since the possible divergences of 
the of the integrand (which occur when k — > 1 and/or when K v — > 1) are compensated by the 
fact that when k — > 1 one has that k' — > as well as when K v — » 1 one has that c(A) — > 
(where K v and c(A) are defined below in the main text). 
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where C(2) is the two dimensional contribution and C(3) the three dimensional 
one. A trivial but cumbersome computation gives 

h h 



(3) 



/(") 

A' 

c(A) 
dc(A) 



A' 
A" 



-4a 



d V [A'c(A)r 



o o 



1 ~ A-) 2 (1 - /(„>) 



3/2 



A"c(A) + (A') 2 ^ 



V 1 - 7h 



(17) 



(sin 2 <p) 



,VA(0)-1 



= fyA, A" = ^A', 7i 



A(0) 
1 

" TV' 



/2 = 



A 2 : 



, , . 1 (2- A) TV 2 

( sin2 <P) v^- 1 - ] Z +v ' , = c(A)A' 



i/2 



v-l 



A 2+v 



(sin 2 <p) 



A 2 +" 
v — 2 



(2- A) (A-l) — 



+ (A-l)^ + ( 2 + ^A-D - (2 _ A) 



= -4a 



1 + (l -fc(/?) 2 sin 2 ^)" (1 - fc(/3) 2 sin 2 



(l-fc(/3) 2 sin 2 0)" 1 kk' 
1 + (1 - fc(/3) 2 sin 2 </>) 



+ ■ 



2fcfc'(a(l-fc(/?) 2 sin 2 </>) Q 1 
(a-1) 



+ 



1 - k{(3) 2 sin 2 (j> 
1 + (1 - k{(if sin 2 0)" 



+ 



*1 ^ 
fc + fc' 



A = 1 



\ 



VAW-1 
A(0) 



sin 2 ip, 1 < A < 2 



where A and c(A) only depend on /3 and </>. It is worth to note that the 
inner integrals on ip can be computed in terms of generalized elliptic functions: 
they are formally similar to the term occurring in the two-dimensional model 
provided one substitute 

k 2 -> K v . 
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Thus, by treating sin 2 <p, 1\ and 1% as constant^, one can rewrite Ct$\ as follows 



C (3) ~ /d<^[A'c(A)] 



I 2 EW (K v )-^E^ (K v ) 



A"c(A) + (A 



,,2 dc(A) 
dA 



(18) 



E {n) (K v ) = I dip 



It is worth to note that E {1 \ £ (2) and E {3) diverge when K v approaches to 1. 
However, in the first two terms 

K v —>1=> c (A) — > 

so that the divergences of the elliptic functions is canceled by the zeros of c (A). 
One can see this as follows 



K^2 \K V -1\ 



n — 2, 3, a\ < a 2 = 1 =>■ 



l| a " 



(v^)^-(l + f) y 



(A - 2) , b n > 0. 



v has to be such that the product 
(c(A)fE^(K u ) ~ 



-i + i 



< oo 



is finite for (A — 2) approaching to zero so that v has to fulfil 

v>v*>\. 

This constraint, as it will be pointed out in the next section, can be interpreted 
as a positive self-consistency check of the present scheme. In the last term in 
Eq. (fl8|) there is a logaritmic divergence when K v approaches to 1 

eM(K v ) ~ log|A-2| 

A^2 



4 Being 1 < N < 2 bounded, the I{ are also bounded and therefore they do not affect 
the critical behavior. Moreover, sin 2 !/? appears in the numerator so that, to understand the 
structure of the possible divergences, sin 2 !/? can be treated as a constant. 
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and this is not compensated by c(A). In fact, the outer integral in <j> smooths 
out this divergence: in the potentially dangerous region such a contribution is 



{K v ) 



A^2 



'log | A 



A->2 



A->2 



d(f) log 



1 — ai sin 



o 2 



< oo, 



where (j>* is the angle such that the argument of the logaritm vanishes and a\ 
and 02 are two real constants. It is then apparent that the divergences comes 
from the terms A" and A': such divergences indeed occur for k approaching to 
1. The term which dominates in this limit is the one coming from A" (such a 
term of A" is the one in which there is no factor k') its contribution being 



C, 



(3) 



1 



(1 - k 2 sin 2 



Therefore, one can choose the parameter a in order to obtain the correct critical 
exponent. It is interesting that the qualitative behaviors for very high and very 
low temperatures are in qualitative agreement with what one expects since the 
novel contribution to the free energy is bounded away from the critical region in 
such a way that the qualitative behavior is not very different from the 2D Ising 
model. The others free parameters (that is, A, v and the overall scale of the 
temperature which has not been explicitly written) do not influence the critical 
region. The overall scale of the temperature should be fixed in such a way that 
k is equal to one at the critical temperature (which, of course, depends on the 
dimension). In 3D it is known (see [TT]) but in any case the present scheme 
(provided the parameters are suitably chosen) seems to be able to describe the 
model in higher dimensions also (in case, with further dressings). 



3 Exploring high temperatures 

It is an interesting results that a very natural ansatze like the one proposed in 
the present approach is able to describe the critical region of the 3D Ising model 
with a single free parameter. This result supports the use of Regge theory (and, 
more generally, of the " old-fashioned" tools of strong interactions based on the 
S matrix and its analytical properties; see, for two detailed review, [4], [TT]) in 
statistical mechanics. The scheme proposed here has a methodological interest 
in itself and could be useful in other situations. For instance, the Sherrington- 
Kirkpatrick solution of the mean field spin glass model (as well as, in many 
respects, the full Parisi solutions) can be naturally interpreted as "dressing" of 
the trivial mean field spin glass model with a highly non trivial expression for 
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the dressing functions giving the spin glass quantities when acting on the mean 
field quantities (see, for instance, [7], [8] and [3]). There are many quantities of 
interest in the mean field spin glass model (such as the magnetization, the free 
energy and so on) which can be written schematically as follows 

m S G(0) = J Dzp SG (z,P)m mf (z) 

where msciP) is the magnetization of the spin glass mean field model, m m f(z) 
is the magnetization of the standard mean field model without disorder and 
p SG {z,(3) is a suitable highly non trivial! dressing function whose explicit ex- 
pression is known (see, for instance, [8]). The above dressing function encodes 
non trivial information about the physical effects of "glassy" disorder. The 
present method suggests a new way to write down suitable ansatze to describe 
non mean field spin glass systems. One could dress the free energy or the magne- 
tization of the 2D Ising model with function(al)s inspired by the corresponding 
dressing function p SG (z,f3) of the mean field model in which some parameters 
should be introduced in order to take into account that one does not expect that 
exactly the same function describes both the mean field and the non mean field 
case. In the spirit of the Regge approach to particles physics, the few param- 
eters which one introduces could be fitted on experimental and/or numerical 
data to provide with a detailed description of the much harder non mean field 
spin glass model. To be more concrete, without entering into the subtle details 
of the theory of mean field spin glass theory (which can be found, for instance, 
in [7], [5] and [3]), the free energy Fg G °f the mean field Sherrington- Kirkpatrick 
model is 

F SG ((3,h)=iu£(log2 + f(0,h;x,P)-^- qx(q)dq) (19) 
x 2 J 

where h is the magnetic field, x is the spin glass functional order parameter 
(which is a non decreasing function of an auxiliary variable q such that both x 
and q take value into [0, 1]). The function f(q, y; x, (3) has a very important role 
and its dependence on the variables q and y is determined by the equation 

d q f+\dlf+ X -^{d v ff=Q (20) 

with final condition 

/(I) J/) = logcosh(/3y). 

The dependence of / on (3 and x is not important for the present discussion. 
In any case, the spin glass formalism gives rise to a non trivial "dressing" of 
the trivial thermodynamic function logcosh(/3y) through the above non linear 
anti-parabolic equation. Therefore, as a first attempt to attach the problem of 

5 To get an idea of the difficulty of the problem, one can think that the discovery of Parisi 
has been done in | 10l . Only recently 5 6 15 , 23 years later, it has been possible to prove 
rigorously that the Parisi solution is the right one to describe the mean field glassy phase. 
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non mean field spin glass systems (and, in particular, the 2 dimensional Ising 
spin glass) one could change the final condition of Eq. ([20)1 and using as final 
condition, instead of the trivial thermodynamic function log cosh(/3y), a function 
similar to the free energy of the 2 dimensional Ising model such as 



f a ,x(l,y) = logcosh(/3y) 



27 



G?X log 



k((3y) 2 sin 2 a;)' 



(in which two parameters a and A have been introduced in order to allow the 
possibility to use both a and A as variational parameters). Indeed, more general 
ansatze are also possible, the basic physical idea remaining the same: according 
to the interpretation of the glassy disorder as a " dressing" of a non glassy system, 
one can change the initial data of Eq. (|20p in such a way to obtain an effective 
description of non mean field glassy systems. I hope to return on this problem 
in a future publication. 

The skeptical reader could simply tells that it is easy to describe the critical 
behavior of a system by writing down a function with some free parameters. On 
the other hand, it is also interesting to notice that the present scheme could also 
be suitable to describe the model in a wider range with the same few parameters. 
Let us consider for instance the small (3 region. The first coefficients of the Taylor 
expansion of for small (3 read 



^ r 3D 



dtfidip 



(T \ 2 M T ^ -L ^ 



^dA f ) 



'(3) 



A" 



<9A : 



(A') 



/(2) ' dA f J v/i -/W 



A'' 



dA 2 



(Af 



3 • 



4ir 2 F, 



2e >,A=1)' 



3D 



A 



/ 00013 

2N 2 



2N 2 



N 3 



1(1)9/3(1(2)) , 1 



2N 



N 3 



(1(3)) 



2N 2 



(1(1)) (1(2) +1(3)) 1 v 



1(1)90 (1(3)) dp 



N 3 



(3) J 



2N 2 



The temperature enters in the above expressions only through the function 
A which is, of course, an even function of (3 through the function k 2 (/3) which 
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vanishes for (3 — together with its odd derivatives. Therefore, in the expansion 
for small (3 for the free energy 



4 2 F (^=i) 

r 3D 



^ a 2 „/3 




fl2n " (2n)!^ 2 " I * ' (21) 



/3=0 



(in which the above coefficient have to be evaluated in j3 = 0) many terms do 
not contribute (such as the one in which at least one factor In\ appear). In 
particular, the term which has the biggest numerical coefficient is the one with 
the highest number of derivatives of Z( 2 ) ■ Generically, one gets 

Such a factor is of order (2n)! which compensates the denominators in Eq. ([2"T]) . 
There are further non vanishing terms in the derivatives of the free energy: the 
number of these further terms increases with n. However, the rate of increasing 
of the coefficients am should increase very slowly with n. In other words, one 
should expect that the ratio 

^±i~in c , oO, t>0, c«l (22) 

«2n 

where c is a very small positive number. The reason is that all the terms in which 
one takes derivative oil/N are zero for (3 = 0. All the terms in which Jm and/or 
7(3) appear without derivatives vanishes for (3 = and so on. Consequently, the 
a,2n cannot increase fast with n. Let us consider the following function 

log(ff(/3)) 

where g is an even function of [3. It is possible to convince oneself that if one 
compute 

c 2n =^"log( 5 (/3))|^ =0 

only in one term it will be present the numerical factor (2n)! compensating 
the denominator in Eq. (|2ip . Moreover, at order 2n, there are at most n non 
vanishing terms contributing to C2„ each of which has a combinatorial factojfl 
smaller than (2n)\. For these reasons, the estimate (|2"21 could hold with t of 
order 10 and c very close to zero. As far as the present scheme is concerned, the 
parameter which is important in this regime is v which is related to how fast 
the free energy changes for small (3. 

It is worth to note that a too large v (such as v > 10) would lead to a fast 
increase of the coefficient am in contrast with the available numerical results 



6 The subleading combinatorial factors are smaller of at least a factor l/2n than the dom- 
inating combinatorial factor. 
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(see, for instance, [T]). As it has been explained in the previous section, a value 
of v lying, for instance, in the interval 

1 < v < 2 

is compatible with the request to have the correct behavior in the critical region. 
Thus, the constraints coming from the critical region are quite consistent with 
the ones coming from the high temperatures region: something which, a priori, 
is far from being obvious. 

These considerations seem to be in qualitative agreement with the observed 
increase of the coefficients computed with numerical methods (see, for instance, 
[1 ). In this scheme the constant A (which has no role near the critical tem- 
perature) measure the " relative importance" of the coefficients of the high tem- 
perature expansion of the 2D Ising model (which is present in the complete 
expression Eq. (fT4|) ) and the genuine 3D term giving rise to the correct critical 
behavior. A also could be extremely important in a numerical fit with the avail- 
able data at high temperature. Unfortunately, the present author is not expert 
in numerical analysis. However, the above considerations suggest that it could 
be possible to choose the parameters A and v in order to achieve a reasonable 
matching with the high temperature numerical results. 

4 Conclusion and perspective 

In this paper a Kallen-Lehman " phcnomenological" approach to 3D Ising model 
has been proposed. The approach is phenomenological in the sense that leads to 
an analytic ansatze for the free energy in which there are free parameters to be 
matched with experimental data and/or to be derived theoretically (as it was 
the case in the Regge-S-matrix approach to strong interactions) . Such a scheme 
seems to be suitable to describe the critical phase. Also the high temperature 
behavior could be captured by the present scheme. It could be also interesting to 
apply the present scheme to non mean field spin glass systems. Unfortunately, 
the present author is not expert in numerical analysis. A careful numerical 
analysis, in case introducing few more parameters, could shed light on the range 
of validity of the present scheme. In any case, the present proposal could also 
have a methodological interest in that similar techniques can be applied in other 
situations. 
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